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We propose a novel mechanism for DCC formation in a first-order chi- 
ral phase transition. In this case the effective potential for the chiral order 
parameter has a local minimum at $ ~ in which the chiral field can be 
"trapped". If the expansion is sufficiently fast a bubble of disoriented chiral 
field can emerge and decouple from the rest of the fireball. The bubble may 
overshoot the mixed phase and subsequently supercool until the barrier disap- 
pears, when the potential resembles that at T = 0. This situation corresponds 
to the initial condition realized in a "quench". Thus, the subsequent align- 
ment in the vacuum direction leads to strong amplification of low momentum 
modes of the pion field. We propose that these DCCs could accompany the 
previously suggested baryon rapidity fluctuations. 

Relativistic heavy-ion collisions might offer the interesting opportunity to study chiral 
symmetry restoration at non-zero temperature and density, which could possibly lead to the 
formation of domains of disoriented chiral condensate (DCC) The strongest amplifi- 

cation of the pion field is obtained for the so-called "quenched" initial condition [^. It is 
assumed that the heatbath is removed instantaneously after restoration of chiral symmetry. 

However, dynamical simulations 1^,^ show that the "quench" does not emerge naturally 
in a heavy ion collision, if the chiral phase transition is second-order or a smooth crossover. 
In this letter, we instead propose a new approach to obtain the "quenched" initial conditions 
naturally in the presence of a first-order phase transition. 

It has been argued that the phase transition for two massless quarks at baryon- 
chemical potential = is second-order which then becomes a smooth crossover for small 
quark masses. On the other hand, a first-order phase transition is predicted for small 
temperatures and large ^. If, indeed, there is a smooth crossover for fi = and non-zero T, 
and a first-order transition for small T and non-zero fi, then the first-order phase transition 
line in the (/U, T) plane must end in a second-order critical point. This point is predicted 
to be at T ~ 100 MeV and ^ ~ 600 MeV. However, some lattice QCD results indicate a 
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first-order transition even at vanishing baryon-chemical potential |^. 

Such temperatures and baryon-chemical potentials can be reached in the central region 
of heavy-ion collisions in the forthcoming Pb(40 AGeV)+Pb experiments at the CERN- 
SPS p], and in the fragmentation regions of more energetic collisions at the CERN-SPS, 
BNL-RHIC, and CERN-LHC (^i ~ 20, 200, 5000 AGeV) 0. Furthermore, fluctuations 
in individual events can also provide rapidity bins with significantly higher /i and lower T 
than on average [1T0|-|13[]. In any case, the dynamical scenario for DCC formation described 



in this Letter applies to the case of a first-order chiral phase transition, and is qualitatively 
independent of the value of /i. Our calculation described below has been performed at /i = 0, 
and the parameters of the Lagrangian (in particular the coupling of the chiral field to the 
heat-bath) have been chosen such as to yield a first-order phase transition. It can be viewed 
as a representative example to illustrate the idea; the basic mechanism works equally well 
also at non-zero yU. 

In case of a first-order transition, the thermodynamical potential as function of the order 
parameter $ exhibits a local minimum at $ = both in the chirally restored phase as well 
as in the mixed phase In high-energy heavy-ion collisions the expansion rate of 



the locally comoving three- volume element can become large This opens the possibility 



that the system can break up into smaller droplets |T2| , [T6| , which might not be able to follow 



an adiabatic expansion. Instead, a bubble can "overshoot" the phase boundary |TT| into the 



low-density broken phase. The chiral field in the bubble is coherent if the bubble-radius is 
on the order of the coherence length or smaller. Close to a first-order phase boundary, the 
chiral field is light and this coherence length is expected to be large. 

Suppose the bubble is created in the restored phase, close to the first-order phase bound- 



ary. Following [jrT| we suppose that the chiral field within the bubble is trapped in the $ ~ 
local minimum of the potential (chiral symmetry is still restored but the field oscillates in 
the false direction), while in previous work |l2| , p!7| it was assumed that chiral symmetry 



breaking had already occured in the bubble. 

The preceding expansion will lead to a velocity profile in the bubble. In other words, the 



bubble will exhibit a Hubble- like expansion with a very large expansion rate [^,0, and 
supercool. During the period of supercooling the chiral field oscillates coherently within the 
whole bubble, while its energy dissipates partly due to friction (coupling to the heat-bath 
p!8|). The local minimum in the effective potential persists until the droplet reaches the 
spinodal line where the potential is close to that at T = (one single minimum). The 
moment where the coherent field "leaps" over the barrier depends on the barrier height and 
the fiuctuation. In a quasi-static situation the field would tunnel to the global minimum. 
However, in a high energy heavy ion collision the local expansion rate is so large (roughly 
10^° times larger than that of the universe at spontaneous chiral symmetry breaking |1T5| ) 
that it is reasonable to assume that tunneling has no time to occur. 
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FIG. 1. Finite temperature effective potential (i.e. the grand canonical potential) within the 
linear cr-model at T = 130 MeV (left) and T = 100 MeV (right). The black dot depicts the average 
chiral field at the corresponding time. 



A second possibility for the chiral field to be kicked over the barrier to the global minimum 
at $ 7^ is via a thermal fiuctuation. Estimates within homogeneous nucleation theory fig 



show that the time-scale for nucleation of critically-sized bubbles with $ 7^ is about the 
same as that needed to reach the spinodal region (^ 2 fm/c in our calculation below). Finite- 
size effects delay bubble nucleation even further |^ : if our entire decoupled bubble has the 
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size of the critically-sized bubble in homogeneous nucleation theory, it can not convert to 
the broken phase via homogeneous nucleation. Instead, it must supercool until it reaches 



the spinodal instability (when the potential resembles that at T = 0), cf. also |jTl|. Thus, 
the most favorable scenario for amplification of the low-momentum modes of the pion field, 
i.e. the "quenched" inital condition, is automatically realized in a natural way if the chiral 
phase transition is first-order. 

Similar to previous studies of DCC formation our scenario also requires a rapid evolution 
out of equilibrium. However, the required 10%-20% supercooling appear much more mod- 
erate than the instantaneous removal of the heat-bath at T ~ Tc, as is necessary to obtain 
the "quenched" initial conditions in a smooth crossover P]. 

After the true vacuum is reached, the coherent chiral field will eventually decay into 



pions due to residual interactions @J2T|. If, at this stage, the heat-bath is still very hot 
and dense, scattering of the DCC-pions with particles from the heat-bath will randomize 



the isospin orientation and spread the momentum distribution |^2[. However, in case of a 
first-order transition with supercooling, the DCC decays at a much lower temperature and 
density of the heat-bath. Therefore, it might be more feasible to detect the DCC-pions. 
These pions will be blue-shifted, though, according to the velocity of the bubble from which 
they emerged. 

Above we discussed the mechanism for DCC formation within a decoupled bubble. In 
principle, however, the same idea can be applied to the entire fireball, thus assuming that 
it supercools and reaches the spinodal instability as a whole. Nevertheless, we have chosen 
to describe how our picture works in a smaller droplet. 

To illustrate the above idea, we applied the linear cr-model coupled to a heat-bath 0,|16 . 
The Lagrangian of the linear sigma model with quark degrees of freedom reads 

£ = q[ird^ - g{a + zj,r ■ n)]q + ^{d^ad^a + d^nd^^) - U{a, vf), (1) 

where the zero temperature potential is 

U{a,^r) = ^ia^ + 7r^-vy-Ha. (2) 



Here q is the light quark field q = {u, d). The scalar field a and the pion field vf = (vri, 7r2, tts) 
together form a chiral field $ = (cr, tt). This Lagrangian is invariant under chiral SUl{2) ® 
SUr{2) transformations if the explicit symmetry breaking term Ha is zero. The parameters 
of the Lagrangian are usually chosen such that the chiral symmetry is spontaneously broken 
in the vacuum and the expectation values of the meson fields are (cr) = and (vf) = 0, 
where /tt = 93 MeV is the pion decay constant. The constant H is fixed by the PC AC 
relation which gives H = fnm'^, where m^r = 138 MeV is the pion mass. Then one finds 
= fl-'^. The sigma mass, = 2X^ fl + m^, which we set to 600 MeV yields ^ 20. 
The coupling to the heat-bath, g = 5.5, is chosen such as to obtain the potential shown in 
Fig. It corresponds to a first-order chiral phase transition. Note that this large value for 
g results in a constituent quark mass at T = of nig = gf^^ ~ 512 MeV. Thus, the nucleon 
mass is too large. However, this is not relevant for our present considerations. A stronger 
coupling leads to a stronger first-order phase transition with a more pronounced barrier and 
even larger constituent quark mass; on the other hand, g = 3.3 fits the nucleon mass in the 
vacuum but results in a smooth crossover 0. 

The Euler-Lagrange equations of motion for the fields, 

d^d>^n + \^[a^ + 7f2 - v^]n = -gpp,, (3) 



are solved self-consistently through the effective quark and antiquark mass, nig = gVa'^ + vr^, 
with the continuity equation for the energy-momentum tensor of the heat-bath, which is 
constituted by the quarks: 

d.T^'' + p^d^nig = 0, (4) 

where ps and ppg are the scalar and pseudoscalar densities, respectively. To solve eqs. (0) we 
employ a second-order leap-frog algorithm, while eqs. (HI) are solved with the RHLLE algo- 



rithm PB[ , assuming that T^'^ is that of an ideal fluid in local thermodynamical equilibrium. 



For more details please refer to 




FIG. 2. The ttq and a field strengths within the forward hght-cone for g = 5.5 and initial condi- 
tion $ = (0.1, 0.07, 0, 0)/7r, d^/dt = 0. The initial expansion scalar was chosen as d-u ~ l/(2fm/c). 

Let us consider the evolution starting close to the minimum of the potential shown to the 
left in Fig. |l]. The evolution of the fields is shown in Fig. |[ The chiral field is "trapped" for 
t ~ 2 fm/c in the local minimum until the temperature drops to the value corresponding to 
the potential to the right (T ~ 100 MeV). At this point the bubble is supercooled by about 
15% and the barrier to the global minimum has almost disappeared. Here, the field starts 
to "accelerate" and rolls towards the true vacuum. Note that the field gets an additional 
"kick" because the bubble starts to collapse under the vacuum pressure which at this 
point, after the system has crossed the first order phase transition line, exceeds that of the 
classical field inside the bubble. For this scenario, we find that the pion field is amplified 
by more than a factor of 100, corresponding to about 18 tto's and 3 cr's for an initial radius 
of 4 fm, see Fig. ^. For the simple quench scenario, where the fields evolve in the zero 
temperature potential U and without coupling to any heat-bath, we get similar results: 21 
TTo's and 4 a's. We computed these numbers as described in ||2^j5|]. 
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FIG. 3. Number of a and ttq mesons produced from the decay of the classical chiral field. The 
solid lines correspond to the self-consistent calculation with first order phase transition (initial con- 
ditions as in Fig. |2|); dotted lines are for the "quenched" initial conditions $ = (0.45, 0.07, 0, 0)/,r, 
d^/dt = 0, y = [/ (this corresponds to roughly the same potential energy at the time when 
"roll-down" begins). 



In summary, we discussed a novel mechanism for DCC formation in a first-order chiral 
phase transition. We study a bubble of disoriented chiral field that decouples from the rest 
of the system before reaching the phase boundary. The bubble supercools and grows due 
to the large volume-expansion rate established during preexpansion (before decoupling of 
the bubble). As the effective potential approaches the T = form, the local minimum at 
$ = disappears. This leads to the "quenched" initial condition. The subsequent alignment 
in the vacuum direction leads to very strong amplification of low momentum modes of the 
pion field. Numerical simulations within the linear cr-model coupled to quarks at non-zero 
baryon-chemical potential are in progress. However, the above general discussion does not 
depend on a specific model for the chiral dynamics. 

If the chiral phase transition is first-order, as is particularly likely the case for finite 
baryon density and if the system breaks up into smaller droplets or bubbles, rapidity 
fluctuations (e.g. of baryon number) can occur, as proposed in |lTTl- [l3i|l6| . If a DCC is 
formed, in coincidence "pion-spikes" of given bubble-isospin could appear in the same pt 
and rapidity range. 
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